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Abstract 



A measuring unit (114) establishes the instantaneous location and movement parameters of a flying object. 
A first extrapolation unit (118) predicts future values of the parameters, if necessary under the assumption 
that one of the parameters stays constant. 

A second extrapolation unit (120,122,124.126.130) considers an extrapolation error as a difference between 
the parameters calculated for a certain time and the parameters received from the object itself for this time, 
minimising the error so that It matches an error model in a continuous run of values. 
USE/ADVANTAGE - Improved accuracy of prediction. Suitable for military use. 
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Beschreibung 

Die Erfindung betriffi eine Vorrichtung zur Oberwa- 
chung von Flugobjekten nach dem Oberbegriff des An- 
spruchl. 

Bei der Oberwachung von Flugobjekten werden hau- 
fig Voraussagen iiber die zukiinftige Flugbahn des Flug- 
objekts benotigt, beispielsweise bei der Luftraumuber- 
wachung, urn abschatzen zu konnen, ob es zu einer Kol- 
lision von Flugobjekten konimt» oder bei der Bekamp- 
fung von Flugobjekten aufgrund der endlichen Flugzeit 
der Bekampfungsprojektiie oder dgl. Hierbei ist es be- 
retts bekannt, eine Extrapolationseinheit einzusetzen, 
welche eine Voraussage uber zukunftige Ortsparameter 
des betreffenden Flugobjekts liefert, zumeist unter der 
Annahme, daB einer der Bewegungsparameter des 
Flugobjekts, insbesondere dessen Beschleunigung, kon- 
stani bleibt. Eine Vorricht gemaB dem Oberbegriff des 
Anspruchs 1, bei der der Extrapolationsfehler durch ma- 
nuell eingegebene zusatzUche Parameter bezuglich des 
Flugobjekts verkleinert wird, ist der US-PS 41 46 780 zu 
entnehmen. 

Die Aufgabe der Erfindung liegt darin, die Prognose- 
genauigkeit fOr den zukunftigen Ort des Flugobjektes 
mil baulich einfachen Mittein zu verbessern. 

Diese Aufgabe wird durch die Merkmale des Kenn- 
zeichens des Anspruchs 1 in Verbindung mit den Merk- 
malen des Oberbegriffs gelost Die erste Extrapolations- 
einheit liefert Prognosewerte, welche ausschlieBlich von 
den gemessenen Ortskoordinaten und der ggf. gemesse- 
nen Geschwindigkeit des Flugobjekts zu einem be- 
stimmten Zeitpunkt abhSlngen. Die zweite Extrapola- 
tionseinheit beriicksichtigt den Extrapolationsfehler, al- 
so die Differenz zwischen den fiir einen Zeitpunkt vor- 
ausberechneten Ortsparametern und den zu diesem 
Zeitpunkt tatsachlich vom Flugobjekt eingenommenen 
Ortsparametern und versucht, diesen Extrapolations- 
fehler zu minimierea Dies geschieht durch fortwahren- 
de Anpassung eines dem zeitlichen Verlauf aufeinander< 
folgender Extrapolationsfehler zugrunde gelegten Feh- 
lermodells an die jeweils auftretenden Extrapolations- 
fehler. 

Der Ausbau der bekannten Vorrichtung in eine erfin- 
dungsgemaB arbeitende Vorrichtung erfordert relativ 
geringen Aufwand. da die bekannte erste Extrapola- 
tionseinheit beibehalten werden kann und lediglich die 
zweite hinzugefugt werden muB. 

Um in einfacher Weise die jeweiligen Extrapolations- 
fehler erhalten zu konnen. werden die MaBnahmen ge- 
maB Anspruch 2 vorgeschlagen. 

Der Bauaufwand verringert sich weiterhin, wenn die 
Extrapolationseinheit gemaB Anspruch 3 einen Pradik- 
tionsfilter (Signalextrapolator) umfaOt anstelle bei- 
spielsweise eines Extrapolationsrechners. 

Als Fehlermodell kommt beispielsweise ein PoJynom 
zweiten oder hoheren Grades infrage. Bevorzugt ist je- 
doch,das Fehlermodell gemaB Anspruch 4. welches ei- 
ner zeitinvarianten iinearen Differentialgleichung zwei- 
ter Ordnung entspricht und bei relativ einfacher techni- 
scherReaiisierungbesondersgule Ergebnisse liefert. 

Eine weitere Verbesserung der Voraussagegenauig- 
keit wird durch die Merkmale des Anspruchs 5 erreicht. 
Aufgrund des relativ kurzen Zeitintervalls. in welchem 
zuruckliegende Extrapolationsfehler bei der Fehlermo- 
dellanpassung beriicksichtigt werden. wird ausgeschlos- 
sen. daB bereits abgeschlossene Flugmandver die Extra- 
polation beeintrachiigen. 

Die MaBnahme gemaB Anspruch 6 fuhri ebenfalls zu 
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ein^r Verbesserung der Voraussagegenauigkeit, da die 
groBeren Absolutbetrage des jeweils ermittelten Extra- 
polationsfehlers bessere Aussagefahigkeit besitzen als 
der in der Nahe des Rauschens liegende Extrapolations- 
5 fehler bei kleinerem Voraussage-Zeitintervall. 

Die Erfindung wird im folgenden an bevorzugten 
Ausfuhrungsbeispielen an Hand der Zeichnung erlau- 
tertEszeigt: 

Fig. 1 eine blockschemaartige, wesentlich vereihfach- 

10 te Darstellung der erfindungsgem&fien Vorrichtung zur 
Oberwachung von Flugobjekten; 

Fig. 2 ein vereinfachtes Blockschaltbild einer in Ver- 
bindung mit einem Feuerleitrechner eingesetzten Aus- 
fuhrungsform der erfindimgsgemaBen Uberwachungs- 

15 vorrichtung und 

Fig. 3 ein Funktionsdiagramm fiir das Blockschaltbild 
gemaB Fig. 2. 

Mit Hilfe der allgemein mit 10 bezeichneten Vorrich- 
tung zur Oberwachung von Flugobjekten gemaB Fig. 1 

20 soil eine mdglichst gute Abschatzung fur einen zukiihf- 
tig von einem Flugobjekt 12 eingenommenen Ort er- 
mdglicht werden. In den Fig. ist mit Indexm angedeutet, 
daB der indizierte Parameterwert aus einer Messung 
resultiert und einen dem MeBzeitpunkt zugeordneten 

25 momentanen Wert angibt und mit Indexb die Tatsache, 
daB der indizierte Parameterwert aufgrund einer Extra- 
polationsberechnung ermittelt worden ist Man erkennt 
in Fig. 1 eine MeO- (und Rechen)einheit 14 zur Messung 
des Ones rm des Flugobjekts 12, beispielsweise mittels 

30 Radar (Radarschirm 15). und zur Abschatzung der mo- 
mentanen Geschwindigkeit Vm sowie der momentanen 
Beschleunigung bm- 

Die erhaltenen Werte werden in einer mit 16 bezeich- 
neten ersten Extrapolationseinheit zugefuhrt, welche 

35 beispielsweise unter der Annahme konstanter Beschleu- 
nigung den Ort rb zu einem zukunftigen Zeitpunkt yor- 
ausberechnet. Der Verbesserung dieses vorausberech- 
neten Wertes rb dient eine zweite Extrapolationseinheit 
18, welche einen mit Arb bezeichneten Extrapolations- 

40 fehler berechnet und mit einem Subtraktionsglied 20 
zufiihrt. welches f rb — Arb bildet. Dieser korrigierte 
Wert f kann weiterverarbeitet werden, beispielsweise 
zur Anzeige des zukunftigen Flugobjektortes auf einem 
Bildschirm. Das der Extrapolation in der zweiten Extra- 

45 polationseinheit 18 zugrundeliegende Fehlermodell 
wird an deh zeitlichen Verlauf der tatsachlich aufireten- 
• den Extrapolationsfehler Arm angepaBt. Der tatsachli- 
che Extrapolationsfehler Arm ist die Abweichung des 
vorausberechneten Ortes rb vom dann tatsachlich vom 

50 Flugobjekt 12 eingenommenen Ort rm- 

In den Fig. 2 und 3 ist ein spezieller Anwendungsfall 
der erfindungsgemaBen Vorrichtung dargestellt, nam- 
lich in Verbindung mit einem Feuerleitrechner einer 
Flugabwehreinheit. Die erfindungsgemaBe Vorrichtung 

55 ist allgemein mit 110 bezeichnet und ist zum einen an 
eine Radarmefieinheit 114 mit Radarschirm 115 ange- 
schlossen sowie uber eine Ausgabeleitung 116 an eine 
nicht dargestellte Flugabwehrbatterie. Die RadarmeB- 
einheit 114 liefert MeBwerte rm an einen Feuerleitrechr 

60 ner 118. welcher die Beschleunigung bm und Richtwin- 
kel ttm (Seite) und pm (Hohe) berechnet. Der Feuerleit- 
rechner 118 ist sowohl an einen Speicher 120 ange- 
schlossen als.auch an ein erstes und ein zweites Subtrak- 
tionsglied 122 und 124. Das erste Subtraktionsglied 122 

€5 ist eingangsseitig mit dem Speicher 120 verbunden und 
ausgangsseitig mit dem Eingang eines Pradiktionsfiliers 
(Signalextrapolator) 126, Der Ausgang des Pradiktions- 
filters 126 wiederum ist mit dem zweiten Subtraktions-. 
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glied 124 verbunden, dessen Ausgang an die erwahnte 
Ausgangsleiiung 116 der Vorrichtung 110 angeschlos- 
sen ist Eine Systemuhr 130 steuert sowohl die Radar- 
meCeinheit 116 wie auch den Feuerleitrechner 118 so- 
wie den Speicher 120. 

Der Funktionsablauf der Anordnung gemaQ Fig. 2 
wird im folgenden an Hand von Fig. 3 erlautert. 

Zu einem willkiirlich herausgegriffenen Zeitpunkt, 
beispielsweise dem Zeitpunkt k— 1 wird eine Radarmes- 
sung durchgefuhrt und der momentane Ort rm sowie die 
momentane Geschwindigkeit Vm und Beschleunigung 
bm (im Feuerleitrechner 1 18) ermitielt. Dies ist im strich- 
lierl angedeuteten Block 132 in Fig. 3 angedeutet. An- 
schlieBend wird durch den Feuerleitrechner 118 ggf. un- 
ter der Annahme gleichbleibender Beschleunigung bm 
der Ort rb des Rugobjekts zur Zeit k' - k + Tk voraus- 
berechnet. Hierbei ist mit k der auf den MeBzeitpunkt 
k— 1 n§chstfolgende MeBzeitpunkt bezeichnet und mit 
Tk die Flugzeit des zur Bekampfung des Flugobjekts 
bestimmten Projektils bis zum Errelchen des Orts rb. In 
einem mit 134 bezeichneten Block des Zeitdiagramms 
gemaB Fig. 3 mit senkrecht nach unten verlaufender 
Zeitachse wird nach Berechnung von rb die Berechnung 
von Richtwinkeln ab und pb (Seite und H6he) vorge- 
nommen. Der vorausberechnete On rb entspricht je- 
doch nicht dem dann tatsachlich erreichten Ort rm. was 
vor allem auf MeBfehler der RadarmeBeinheit 114 zu- 
ruckzufuhren ist und ggf. auch auf plotzliche Bewe- 
gungsanderung des Flugobjekts 12, Es hat sich nun ge- 
zeigt, daB in vielen Fallen diese Extrapolationsfehler 
einen charakteristischen zeitlichen Verlauf aufweisen. 
welcher in befriedigender Weise durch eine stochastisch 
erregte ged^mpfte Schwingung angenahert werden 
kann mit den Parametern Eigenfrequenz, Dampfung 
und Starke der Erregung. An Hand der jeweils gemesse- 
nen Extrapolationsfehler in einem Zeitintervall, welcher 
vom momentanen Zeitpunkt jeweils einige Sekunden 
zuriickreicht, lassen sich die angegebenen Modellpara- 
meter in ausreichender Genauigkeit ermitteln. An Hand 
des so angeps^Bten Modells lassen sich die f iir den Zeit- 
punkt t = k zu erwartenden Extrapolationsfehler vor- 
ausberechnen, welche im vorliegenden Falle Rich twin- 
kel-Korrektiiren Aab und Apb bilden. Diese im Block 
136 gemaB Fig. 3 bzw. im Pradiktionsfilter 1* gemaB 
Fig. 2 gebildelen Korrekturwerte werden im Subtrak- 
tionsglied 124 bzw. im Block 138 von den vorausberech- 
neten Werten ab und pb abgezogen, worauf man die 
korrigierten Winkelwerte a und p erhalt. Diese korri- 
gierten Werte werden iiber die Leitung 1 16 der Flugab- 
wehrbatterie zugefuhrt. 

Im Block 140 gemaB Fig. 3 wird nun mit Hilfe einer 
entsprechenden Radarmessung zum Zeitpunkt k' der 
vom Flugobjeki zum angenommenen Trefferzeitpunkt 
k' = k + Tk tatsachlich eingenommene Ort rm festge- 
stellt und hieraus mittels des Feuerleitrechners 1 18 eine 
Berechnung der Richtwinkel am und Pm vorgenommen, 
d h. derjenige Richtwinkel, welche zu einem Treffer am 
Ort rm fuhren wiirden. 

Als nachstes wird im Block 142 gemaB Fig. 3 mit Hilfe 
des Subtraktionsglieds 122 in Fig. 2 der tatsachlich auf- 
tretende Extrapolationsfehler Aam = Om — ab und 
APm = Pm - 3b berechnet. Es ist klar. daB der Speicher 
120 bis zu diesem Zeitpunkt k' die bereits zum Zeit- 
punkt k ermittelten Werte ab und Pb abgespeicheri hat. 

Diese Extrapolationsfehler Aam und APm werden nun 
zur Beendigung dieser Zeiischleife dem Block 136 in 
Fig. 3 entsprechend dem Pradiktionsfilter 126 in Fig. 2 
zugefuhrt zur Anpassuiig der Modellparameter. Es be- 



glnnt eine neue Schleife mit einer weiteren Radarmes- 
sung im Block 132 gemMO Fig. 3 zu einem nachstfolgen- 
den Zeitpunkt. Dieser Zeitpunkt kann spater liegen als 
der Zeitpunkt k' oder auch friiher (z. B. zum Zeitpunkt 
5 k), wobei dann dementsprechend durch entsprechende 
Ausbildung des Speichers 120 und entsprechende 
gleichzeitige Steuerung durch die Systemuhr 130 dafiir 
gesorgt werden muB» daB dementsprechend viele Werte 
ab und Pb im Speicher 120 zwischengespeichert und 
10 jeweils dann an das Subtraktionsglied 122 abgegeben 
werden. wenn die zugeordneten Winkelwerte ttm und pm 
zur Verfiigung stehen. 

Die Anpassung der Parameter des Fehlermodells an 
die jeweils gemessenen Extrapolationsfehler Aom und 
15 APm sowie die Vorausberechnung zukunftiger Extrapo- 
lationsfehler Aab und APb wird im Pradiktionsfilter 126 
vorgenommen. Dieser kann auch als Signalextrapolator. 
Vorhersagefilter, PrSdiktor oder Datenextrapolator be- 
zeichnet werden. Mit Hilfe von Extrapolationsfiltern. 
20 wie diese bei der Digital-Analog- Wandlung haufig ein- 
gesetzt werden (s. beispielsweise Taschenbuch "Elektro- 
technik" Karl Hanser-Verlag, Munchen — Wien 1978, 
Band 111, ab Seite 898), laBt sich aus der bekannten Ver- 
gangenheit einer Zahlenfolge das zugehorige kontinu- 
25 ierltche analoge zeitabhangige Signal extrapolieren, und 
beispielsweise ein zukQnfliger Funktionswert berech- 
nen. Der. Extrapolation konnen beispielsweise Polyno- 
me zugrunde gelegt werden oder, wie im vorliegenden 
Falle, eine stochastisch erregte gedampfte Schwingung. 
30 Dem Entwurf des Pradiktionsfilters 126 liegt ein Zu- 
standsmodell einer zeitinvarianten, linearen. stochasti- 
schen Differenzialgleichung zweiter Ordnung zugrunde. 
Mit Hilfe des Kalman-Filter-Algorithmus (s. z. B. 
Ch. Landgraf: stochastische lineare Systeme, TU Berlin, 
35 Brennpunkt Kybemetik 1979) laBt sich der Pradiktions- 
filter 126 mit den gewunschten Eigenschaften entwer- 
fen. Am Eingang des Pradiktionsfilters liegen die gemes- 
senen Extrapolationsfehler Aam und APm an; das Aus- 
gangssignal des Pradiktionsfilters bilden die vorausbe- 
40 rechneten Extrapolationsfehler Aab und APb- Die Band- 
breite des Pradiktionsfilters 126 ist derart festgelegt, 
daB EinfluB welter zuruckliegender gemessener Sub- 
traktionsfehier derart rasch abfallt, daB im wesentlichen 
lediglich gemessene Extrapolationsfehler aus einem im 
45 Sekundenbereich liegenden Zeitintervall sich auf die 
Anpassung der Modellparameter (Kreisfrequenz, 
Dampfungskonstante und ErregungsgroBe) auswirken. 

Bei der Ausfuhrungsform gemaB Fig. 2 und 3 werden 
aus den gemessenen und extrapolierten Orten rm bzw. 
50 Tb jeweils Richtwinkel a und p berechnet, die wiederum 
zur Extrapolationsfehlerminimierung in den Blocken 
142 und IMverwendet werden. Es ist offensichtlich, daB 
in anderen Anwendungsfallen, beispielsweise in dem in 
Fig. 1 dargestellten Fall, auf die Zwischenrechnung (Be- 
55 rechnung der Richtwinkel aus den Orten) verzichtet 
werden kann, unter entsprechender Vereinfachung des 
Zeitdiagramms gemaB Fig. 3. 

GemaB den Figuren ist im Vorstehenden jeweils un- 
ter r der Ortssektor F* zu verstehen und dementspre- 
60 chend unter v und b jeweils der Geschwindigkeitsvektor 
V**' und der Beschleunigungsvektor B*. 

Patentanspruche 

65 1. Vorrichtung zur Oberwachung von Flugobjekten 
(12) mit einer MeBeinheit (14; 114) zur Feststellung 
von momentanen Oris-Bewegungsparametern (rm; 
am. Pm) eines Flugobjektes (12) und mit einer ersten 
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Extrapolationseinheit (16; 118) zur Voraussage zu- 
kunf tiger Ortsparameter (fb; ab, Pb), ggf- unter der 
Annahme, daB einer der Bewegungsparameter (b) 
des Flugobjekts (12) konstant bleibt, dadurch ga- 
le nnzeichnet, daB cine zwette Extrapolationsein- 5 
heit (18; 20; 120, 122, 124, 126, 130) vorgesehen ist, 
welche den Extrapolationsfehler (Ar; Act, AP) als 
Differenz zwischen den fur einen Zeitpunkt (k') 
vorausberechneten Ortsparametern (rb; ab. Pb) und 
den zu diesem Zeitpunkt tatsachlich vom Flugob- to 
jekt (12) eingenommenen Ortsparametern (rm; am, 
Pm) berucksichtigt und diesen Extrapolationsfehler 
(Ar; Aa; Ap) dadurch minimiert, daB es ein dem 
zeitlichen Verlauf aufeinanderfolgender Extrapola- 
tionsfehier zugrunde gelegtes Fehlermodell fort- 15 
laufend anpaBt 

2. Vorrichtung nach Anspruch 1, gekennzeichnet 
durch einen an die erste Extrapolationseinheit (118) 
angeschlossenen Zwischenspeicher (120) fiir die 
vorausberechneten Ortsparameter (ab, Pb) sowie 20 
ein an den Zwischenspeicher (120) angeschlossenes 
Subtraktionsglied (122) zur Ermittlung des Extra- 
polationsfehlers (Aam, APm). . 

3. Vorrichtung nach Anspruch 1 oder 2, dadurch 
gekennzeichnet daB die zweite Extrapolationsein- 25 
heit (120 bis 130) einen dem Fehlermodell entspre- 
chend ausgelegten, ggf. an den Ausgang des Sub- 
traktionsglieds (122) angeschlossenen Pradiktions- 
filter (126) umfaBt 

4. Vorrichtung nach einem der vprhergehenden 30 
Anspruche, dadurch gekennzeichnet, daB als Feh- 
lermodell eine ged^mpfte Schwingung. vorzugs- 
weise mit stochastischer Erregung, eingesetzt ist. 

5. Vorrichtung nach einem der vorhergehenden 
Anspruche, dadurch gekennzeichnet, dafi die zwei- 35 
te Extrapolationseinheit (120 bis 130) zur Optimie- 
rung des Fehlermodells im wesentlichen lediglich 
zuruckliegende gemessene Extrapolationsfehler in 
einem Zeitintervall im Bereich zwischen 0,1 Sek 
und 50 Sek.. vorzugsweise 1 bis 10 Sek. berucksich- 40 
tigU ggf. durch entsprechende Festlegung der 
Bandbreite des Pradiktionsfilters(126). 

6. Vorrichtung nach einem der vorhergehenden 
Anspruche, dadurch gekennzeichnet, daB das Vor- 
aussage-Zeitintervall (k' — k) groQ ist gegenuber 45 
dem Zeitabstand aufeinanderfolgender Ortsmes- 
sungen und vorzugsweise im Bereich zwischen 0,1 
und 20 Sek. am besten zwischen I und 10 Selcliegt. 

7. Vorrichtung nach einem der vorhergehenden 
Anspruche mit einem die erste Extrapolationsein- 50 
heit umfassenden Feuerleitrechner (118) fur die 
Ortsparameter bildende Richtwinkel (ctb und Pb), 
dadurch gekennzeichnet, daB die zweite Extrapola- 
tionseinheit (120 bis 130) die Extrapolationsfehler 
(Aam» APm) als Differenz zwischen den errechneten 55 
Richtwinkeln (ab und pb) fur einen Trefferort zu 
einem bestimmten Treffzeitpunkt und den Richt- 
winkeln (am. Pm), welche dem zum Treffzeitpunkt 
vom Flugobjekt (12) tatsachlich eingenommenen 
Ort zugeordnet sind, minimiert. 60 
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